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1 Introduction 


The standard approach of quantum statistical mechanics uses the density op- 
erator p to describe the (mixed) state of the physical system of interest. Since p is 
an operator in the Hilbert space, we usually need some representation to perform 
the practical calculations. There are many equivalent representations in the lit- 
erature, e.g., the coordinate representation, P-representation, Q-representation, 
Fock space representation, Wigner function and characteristic function (Chi func- 
tion hereafter), etc. In this paper we will use the last two representations since 
they provide a “phase-space picture” for the quantum-mechanical problems [1]. 

In quantum equilibrium statistical mechanics, a system (A) immersed in a heat 
bath (B) with temperature T is described by the canonical ensemble. According 
to ensemble theory, the density operator is: 

_ exp (—/?#) m 

Tr[exp(— /?//)] ’ 

where /? = ^ is the inverse temperature and H is the Hamiltonian of (A). The 
structure of (B) and the interaction between (A) and (B) are irrelevant to the 
density operator. If H is (inhomogeneously) quadratic and with a finite number of 
degrees of freedom, the density operator will correspond to a finite-mode thermal 
Squeezed Coherent State (SqCS) [2]. 

In quantum non-equilibrium statistical mechanics, ensemble theory is no 
longer valid and we have to build a model for the heat bath (B) and consider 
(A)+(B) as a total system. The total Hamiltonian then contains three parts — 
the Hamiltonian of (A) and (B) and the interaction Hamiltonian. 

It is well known that the number of degrees of freedom of a heat bath must 
be infinite (the thermodynamic limit), otherwise, due to the Poincare recurrence 
theorem there will be no phenomena such as approach to equilibrium, damping 
or dissipation. The simplest model of a heat bath is an assembly of infinitely 
many harmonic oscillators with linear couplings to (A). In this kind of model, 
the total Hamiltonian is quadratic if the Hamiltonian of (A) is quadratic. Since 
for quadratic Hamiltonians we have a phase-space picture of quantum mechanics 
with the help of Wigner and/or Chi function, we can construct an infinite-mode 
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(pure) SqCS for the total system in an infinite-dimensional phase space using 
these functions. After reduction, i.e., ignoring the heat bath but keeping its “in- 
fluence”, we will get a finite-mode SqCS for (A). In the limit as time approaches 
infinity, it can be shown that (A) will approach equilibrium, and the finite-mode 
SqCS will become a thermal SqCS consistent with the fluctuation-dissipation 
theorem [3, 4, 5, 6, 7, 8]. 

In this paper we introduce a geometric interpretation of these non-equilibrium 
phenomena via the Chi-function representation of infinite-mode SqCS. In Sec. 2 
notations, conventions and a lemma on matrix are introduced for the mathematics 
used in this paper. In Sec. 3 we study finite-mode SqCS’s by Wigner and Chi 
functions and then extend them to infinite mode. In Sec. 4 we use the quantum 
Brownian motion as an example to illustrate geometrical reduction in phase space. 


2 Mathematical Preliminaries 

Throughout this paper, h is set equal to 1; “T” denotes the transpose of a 
matrix and T”denotes the inverse of the transpose of a matrix. The phys- 
ical system under consideration is of N = n + 1 degrees of freedom, where n 
is either finite or equal to infinity. The symbols x = (x 0 ,^ 1 ,# 2 , . . . , x n ) and 
k = (k°, k l y k 2 , . . . , k n ) denote the TV-dimensional canonical coordinate and mo- 
mentum respectively, thus z = (x; k) is a vector in 2 A— dimensional phase space. 
q and p denote TT-dimensional position and momentum operators corresponding 
to the canonical variables x and k. 

The metaplectic group Mp(2N,R) is an N(2N + l)-dimensional Lie group. 
It is the quantum analogue of symplectic group Sp(2N,R) . The elements of the 
Lie algebra of Mp(2N,R) can be organized as anti-hermitian operators in the 
following form: 

i " 

= jE [ Q oW + PijPiPj + 7 iMiPj + Mi)] 

= ^{q\p)J™{q\p) r , ( 2 ) 
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where otij = atji, 0ij = 0 JX and 




G sp(2N,r) 


is a 2 N x 2 N real matrix [9], while 


J = 




1 = N x N unit matrix. 


( 3 ) 

( 4 ) 


The Lie algebra of Mp(2N,R) is isomorphic to that of Sp(2N,R). The action 
of exp[$(m)] € Mp(2N,R) on ( q\ p) is: 

exp[$(m)](|;p) T exp[-$(m)] = exp(-m)(|;|) T , (5) 


where exp(— m) G Sp(2N,R). 


Lemma [10] 

If M is a symmetric and positive definite 2 N x 2N matrix, then there exist 
two matrices Si,Si G Sp(2N,R) (but not unique), such that 

M = * (o °i) s - < 6 > 

where u> = diag(u>b,t<;i,u> 2 , . . . ,w n ), u>j > 0 for all j, and 



Remarks: 

(1) S G Sp(2N,R) if and only if S T JS = J by definition. 

(2) u>j’s are not eigenvalues of M in general. We will call them the “symplectic 
eigenvalues” of matrix M. 

(3) The eigenvalues of JM are iron’s, hence we can calculate ojj's from JM 
as an ordinary eigenvalue problem. 

(4) If the matrix Cj corresponds to a 2-dimensional rotation on the 
( Xj,kj ) plane, then 




( 8 ) 


Therefore Si in (6) can be replaced by CjS\ and hence is not unique. 
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3 Squeezed Coherent States in Phase Space 

3.1 Wigner Function 

The Wigner function of an N-mode density operator p is defined as [11, 12]: 
W(x\ k) = ir~ N [ dyexp(2ik-y)p(x — y,x + y), (9) 

where p(x, x') is the coordinate representation of the density operator p. 

The Wigner function is real and normalized by definition: 

r dxdkw(x-,k) = i. (io) 

J — OO 

However, it is not always positive— definite and is thus called the quasi— probability 
distribution function over the “phase space” z — (x\k). 

If the density operator is an exponential of a quadratic form of position and 
momentum operators, then the Wigner function is a Gaussian distribution in z: 

W{z) = CVexp[— (5*— z c )M(z - z c ) T ]. (11) 

where Cn = ir~ N iJdet(M) is the normalization constant, z c is a constant vector in 
the 21V— dimensional phase space, and M is a symmetric, positive definite matrix 
with all its symplectic eigenvalues smaller or equal to 1. (Otherwise (11) will not 
correspond to a physical state.) The Gaussian Wigner function (11) corresponds 
to the multimode thermal SqCS in general [2]. It contains the ordinary coherent 
states (when M is a unit matrix) and the ordinary SqCS (when all the symplectic 
eigenvalues of M equal to 1) as special cases. 

The “Wigner ellipsoid” of (11) is defined as: 

(?-z c )M(z-z c ) T = 1, (12) 

which is an ellipsoid in the 2,/V-dimensional phase space with its center at z c and 
its shape determined by M. We can take the Wigner ellipsoid as a geometric 
representation of the Gaussian Wigner function. 
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3.2 Characteristic Function (Chi Function) 

The Chi function of a density operator p is defined as: 

X(x; k) = Tr[pD(-x] -%)], (13) 

It can be shown that Chi function is the symplectic Fourier transformation of 
the Wigner function: 

x(x; k) = f dx'dk'W (x'\ k ') exp[— i{x' • k — k' ■ x)]. (14) 

J — OO 

The normalization condition of the Wigner function corresponds to x(0; 0) = 1 
in the Chi function. Since the operator D{—x\ —k) is unitary, k) is complex 
in general. 

The Chi function which corresponds to the Gaussian Wigner function (11) is 
also Gaussian: 

x(z) = exp[— ^zJM' 1 J T z T + izJzJ]. (15) 

Analogue to the Wigner ellipsoid, we can also define the “Chi ellipsoid” for a 
Gaussian Chi function as: 


(z- z c )JM~ l J r (z — z c ) T = 1. (16) 

The center of Chi ellipsoid is the same as that of the Wigner ellipsoid, while the 
shape of this ellipsoid is determined by the matrix JM~ l J T . 


3.3 Mean Vector and Covariance Matrix 


For an jV-mode (mixed) state with the density operator p, the mean vector 
is defined as ((<?); (p)) in the 2A r -dimensional phase space, where (<?,) = Tr(pq t ), 
etc. The covariance matrix is defined as a 2 N x 2N matrix of the form: 



U* = ((?* - (&»(& - <<7j») = (Wj) ~ <?«)(9j), 
Vij = (( Pi - ( Pi))(pj - ( Pi ))) = (PiPj) - (Pi)(Pj), 


(17) 

(18) 
(19) 
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Qij = !<(& - <«))(&* - <Pi» + (Pi - <p;»(& - <&»> 

= <^(9iPj+Pi9i))-(*)(Pi)- ( 20 ) 

For the Gaussian Wigner function (11) or Gaussian Chi function (15), the mean 
vector is z c , and the covariance matrix takes the form: 



Therefore (15) can be re-written as: 

= - ® T )* T + «■*?]• (22) 

3.4 Time Evolution of Wigner and Chi Functions 

Consider an iV-mode Hamiltonian: 

H = \(hv)K (?;?) T , (23) 

where A" is a 2 N x 2 N positive definite symmetric matrix . According to the 
lemma introduced in Sec. 1, this kind of Hamiltonian can be transformed into 
the following form: 

H = \(IJ)S T J)s(?;p), (24) 

where w = diag(a\),w 1 , . . . ,u> n ) and the Wj’s are symplectic eigenvalues of the 
matrix K. The time-evolution operator exp (—itH) is an element in Mp(2N, 
R) and the time evolution of (q\p) is a special case of (5): 

exp(itH)(q-,p) r exp(-itH) = R(t)(|;p) T , (25) 

where 

„ i ( cosfu >t) uj 1 sin(u>t) 

R(t) = exp (tJK) =S~ l [ : \ ' 

\ — w sin(u >t) cos(u>t) 

is an element in Sp(2N,R). {f?(t)| t €R) forms a one-parameter subgroup of 

Sp(2N,R) describing the phase flow in the 2,/V-dimensional classical phase space: 

(x(t); k(t)) T = R(t)(x( 0); k( 0)) T . (27) 
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It is well known that for the quadratic Hamiltonian (23), the time evolution 
of Wigner function and Chi function follow exactly this phase flow: 


W(z\t) = C'vexpH?- z c (t))R(t)- T MR(t)~'( S - ? C (()) T ], 
X(z\t) = exp[-jM(i)- T (_ V Q _ y T )flW'? T +>?^7(‘)] 

= exp l-2H-0(t) U{t) ) Z +,zJz ‘W' 


( 28 ) 


(29) 


where zj(t ) = R(t)z { 


<• 


3.5 Reduction of Multimode Squeezed Coherent States 

Consider the quantum system (A)+(B) discussed in Sec. 1 whose density 
operator is p A B ■ If we reduce this system by ignoring (B), the expectation value 
of an operator 0 A which corresponds to a measurement on (A) will become: 

{Oa) = Tt\Pa6aI (30) 


where p A = Tr( B )(pAB ) is a well-defined reduced density operator for (A) which 
contains the “influence” of (B) on (A), Tr^j represents the “partial trace” which 
only takes trace with respect to the degrees of freedom of (B). 

If the Wigner function W{x A ,x B \ k A ,k B ) corresponds to the original density 
operator p A B , then the reduced Wigner function corresponding to p A will be [12]: 


W A (x A ;k A ) = f dx B dk B W(x A ,x B ] k A ,k B ). 
J — OO 


(31) 


As for the Chi function, if X(x A ,x B ; k A , k B ) corresponds to p A B, the reduced Chi 
function corresponding to p A will take the form: 


X A {x A ]k A ) = X(x A ,0-,k A ,0), 


(32) 


which is a restriction of the original X(x A \x Bl k A \k B ) to a subspace in the 2N- 
dimensional phase space. From the mathematical point of view, it is easier to use 
the Chi function to perform the reduction. 
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Now let us use TV- mode to one-mode reduction as an example. For a given 
TV-mode Gaussian Chi function (22), we want to make a reduction by ignoring 
all the degrees of freedom which correspond to modes 1,2 , ... ,n and leave only 
the 0-th mode. Without any calculation, we can write down the reduced Chi 
function directly : 

X(i°, k°) = exp[-^(x°, k°) ^ (x°, k°) T + i(x°k° - *qx°)], (33) 

which is a one-mode Gaussian Chi function. 

The geometrical interpretation of this reduction process is cutting the original 
Chi ellipsoid in the 2 TV-dimensional phase space by a “shifted (x°, fc°) plane” — the 
plane which is parallel to (x°, k°) plane and passes through the center of the Chi 
ellipsoid. The section on the Chi ellipsoid gives the “Chi ellipse” on the shifted 
(x°, k°) plane which represents the reduced one-mode Gaussian Chi function. A 
schematic graph of this geometrical reduction is shown in Fig. 1. 

3.6 Infinite-Mode Squeezed Coherent States 

The infinite-mode SqCS is a naive generalization of finite-mode SqCS. Com- 
paring the three equivalent representations of finite-mode SqCS’s: (11) , (15) and 
(22) , we see that (22) can be directly generalized to infinite mode without any am- 
biguity or convergence problem. So we will take (22) in the infinite— dimensional 
phase space as the definition of infinite-mode SqCS, all the above formulas which 
involve (22) can be applied to infinite-mode case. 

4 Quantum Brownian Motion 

In this section we shall study quantum Brownian motion of a harmonic os- 
cillator. The physical picture is a quantum harmonic oscillator immersed in a 
dissipative heat bath. In classical statistical mechanics, this problem can be 
studied via the Langevin equation: 

M X +M7 X +MQ. 2 X = 0, (34) 
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Fig. 1. Reduction as a geometrical operation in phase space. 
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where X, M and fi are the coordinate, mass and frequency of the oscillator, and 
My is the friction constant. 

The quantum analogue of the Langevin equation can be achieved by several 
quantum-mechanical heat-bath models, e.g., the FKM model [3], linear cou- 
pling model [4, 5], independent-oscillator model [13], etc. Actually it can be 
proved that all these models are equivalent [14]. In this paper we will use the 
independent-oscillator model since it is the simplest and most intuitive. 

i 


4.1 Independent— Oscillator Heat— Bath Model 

Consider the Brownian particle to be a harmonic oscillator immersed in a dis- 
sipative heat bath with inverse temperature /?. Using the independent-oscillator 
heat-bath model, the total Hamiltonian of the system is [13]: 


H = 


P 2 
2M 


+ -Mfl 2 Q 2 + + o m * u, « 2 (9i — <5) 2 ]> 


2rrii 2 


(35) 


where Q and P are the position and momentum operators of the Brownian par- 
ticle; qi, pi, ra< and u>, are the position operator, momentum operator, mass and 
frequency of the i— th heat— bath oscillator, i = 1,2,3, .. .. This Hamiltonian is a 
special case of (23). 

It can be proved that in order to make the Brownian particle satisfy the 
quantum Langevin equation: 


M 0 +My £ +M&Q = 0, 


(36) 


the spectral distribution of heat-bath oscillators must obey: 

n 2 

T: mi<J*8(u — u>i) = —My. 

»=i ^ 


(37) 


4.2 Quantum Brown Motion in Phase Space 

In the following we will study time evolution of the Brownian particle by the 
reduced Chi function. The initial condition is chosen to be p = Pa® Pb, where pa 
is the initial density operator of the Brownian particle which corresponds to an 
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arbitrary Gaussian Wigner/Chi function, and ps is the initial density operator 
of heat bath which is in thermal equilibrium at the inverse temperature /?. Since 
the detailed calculations can be obtained by combining the calculations in [6] and 
[13], here we will only discuss the result and the geometrical interpretation. 

Let the degree of freedom of the Brownian particle correspond to the O-th 
mode, and those of the heat bath correspond to other modes. The infinite-mode 
Chi function for the initial condition state is (22) with the following parameters: 
Uoo, Ho, Qoo and 2 °, which correspond to the initial conditions of the Brownian 
particle, are arbitrary; z c has only two non— zero components corresponding to 
2 ® since the mean vectors for all heat— bath oscillators equal to zero; and other 

C 

elements in the covariance matrix are: 


Uij = coth(-/3u>j)£ij, 

(38) 

Vij = ^rrnui coth(^u»i)<5,j, 

(39) 

; = H = Hi = Qio = Qoj = Qij = o, 

(40) 


for all ij =1,2,3, — 

Combining (29) and (33), we get the time-dependent reduced Chi function of 
the Brownian particle (the index 0 for the canonical coordinates is suppressed): 

x(x,M) = exp[-~(x,fc) ^ Uw{t)^) ( x ^ k ) T+i ( xkc ( t )~ kx ^ t ^- ( 41 ) 

It is easy to write down the corresponding Wigner function by comparing (11) 
and (15): 

W(x, k-, t) = C\ exp[— (x - x c (t),k - k c (t))M(t)(x - xjt), k - k c (t)) T ] t (42) 


(43) 


Where ..... 1 ( Ho (() -Ooo(<A 

{ 2[[/„o(<)Voo(<) - Coo(0] V ~Qoo(t) Uoo(t) ) ' 

Unlike ordinary reduction methods [5, 6, 7], we obtained this reduced Wigner 
function without using integration over the heat— bath degrees of freedom. 

Comparing (41) with (42) and (43), we see that at any moment the Wigner 
ellipse and the Chi ellipse are similar and their areas inversely proportional to 
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each other. (Although both areas are time-dependent in general.) When time 
approaches infinity, the Brownian particle will approach the equilibrium state 
which is independent of its initial condition and consistent with the fluctuation- 
dissipation theorem. In Fig. 2, we plot the time evolution of Wigner ellipse and 
Chi ellipse of the Brownian particle in phase space. 


5 Conclusion and Outlook 

The method and result discussed in this paper are valid as long as: (1) The ini- 
tial state of (A) is a finite-mode (not necessary one-mode) Gaussian Wigner/Chi 
function. (2) The Hamiltonian of (A) is quadratic and with finite number of 
degrees of freedom. 

If (1) no longer holds, then we will not be able to use an ellipsoid in phase space 
to represent the state. However, the phase-space picture continues to be valid 
since time evolution of the Wigner/Chi function will still follow the phase flow in 
classical phase space. On the contrary if (2) is not true, e.g., as in quantum tun- 
neling problems, then time evolution of the Wigner/Chi function will not follow 
the phase flow exactly and the phase-space picture will fail. In order to relieve 
this limitation, some authors introduced the idea of “effective potential” [15, 16] 
so that time evolution of the Wigner/Chi function can be still expressed in terms 
of the phase flow. Integration of this modified phase-space picture with the 
dissipation mechanism is an interesting question and worth pursuing further. 
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Fig. 2. Time evolution of Wigner ellipse (solid) and Chi ellipse 
(dotted) of the Brownian particle. Note that the area of the former 
is always larger than that of the latter. 
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